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This theoretical paper is to advance a phenomenological, quasi thermodynamic approach to the dynamics 
of classical liquids which uses the Levy distribution of probability theory. Doubts from the chemical physics 
community about the application of its unusual properties to this field are tried to be removed. In particular, 
to understand the preponderant component of the Levy sum for Glarum Levy defects and Fischer speckles, the 
classical mathematical proof [D. A. Darling, Trans. Amer. Math. Soc. 73, 95 (1952)] for the existence and the 
influence of this component is accompanied by addition of physical arguments related to these defects. It is tried 
to explain an underlying fluctuating spatial pattern of free volume with weak contrast and a pattern of mobility 
with strong contrast, and to explain the characteristic lengths for the main transition and the Fischer modes. The 
structure of the relaxation chart (dynamic glass transition) and several properties of, and relations between, the 
slower dispersion zones therein, are reviewed for classical glassforming liquids of moderate complexity. For the 
main transition, the preponderant component is pushed in the midst of the defect and induces the molecule to 
its diffusion step across the cage door of the next neighbors. An Experimentum Crucis for an indirect proof of 
the existence of defects — via characteristic lengths — is also described. 
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I. INTRODUCTION 

The Levy sum for stable distributions with Levy exponent a < 1 has a preponderant component 1 - 2,3 . An example of its 
physical relevance was recently published for Laser cooling 4 . The reviews for dynamics in classical molecular liquids are 
sometimes accompanied by catalogs of open questions 5 ' 6 - 7 ' 8 . The relevance of Levy distributions for this dynamics 9 may 
be discussed in a thermodynamic frame of fluctuating space-time pattern for free volume, where islands of mobility 10 from 
temporary concentrations of free volume can be formed by using a general length/mobility scaling 11 . Such spots in the pattern 
are called Levy defects, if they can be explained by an appropriate Levy sum, in particular with a preponderant component. A 
well known model is Glarum's defect diffusion model 12 which was later reviewed 13 by means of Levy flights, cf. also RefJ^. 
The more later variant of Bendler, Fontanella, and Shlesinger 15 is aimed to interaction of several or many defects forming 
clusters as basis for the glass transition. On the contrary, our variant 9 is aimed to an explanation of molecular cooperativity of 
the dynamic glass transition from one defect alone. An experimental indication for defects in larger spatial scales are Fischer's 
speckles^. 

If confirmed, the Levy defect pattern would be a surprising phenomenon: Upon some "primary democratic" conditions, Levy 
statistics can promote shaping. 

The basic physical assumption of our approach is that the Levy sum components can be related to spatially separable, 
dynamically independent subsystems or partial systems. (Such a relation is not necessary, which means that there are Levy 
situations where spatial aspects play a minor role. An example is the treatment of Laser cooling in Ref. 4 , where a trap in 
the momentum space is considered, and the Levy distribution is related to lifetimes such as trapping or escaping times.) Our 
assumption is completed by a general space vs. mobility scaling of fluctuation modes to treat the spatial structure of the 
defects. Our approach is supported by a Representativeness Theorem, showing inversely that the distribution of representative 
subsystems corresponds to a Levy sum. An Experimentum Crucis for Levy defects of the main transition in classical liquids is 
suggested that uses the relevant characteristic lengths. 



* * * 

As a rule, experimental retardation in glass forming liquids can well be adjusted by a Kohlrausch functio n 18 ' 19 ' 2Q i 21 in the time 
t domain, 

retardation ~ correlation function ~ exp(— at a ) , a > 0, a = 0.4. ..1.0. (1) 
The proportionality to a correlation function is a consequence of the classical fluctuation dissipation theorem, FDT. 

On the other hand, such "stretched exponentials" are important in probability theory: They are characteristic functions 
(Fourier transforms) for stable (Levy) distributions with exponents a < 1. If the stretched exponentials Eq.(Q} are really 
indications for a Levy distribution, then the latter is some kind of (inverse) Fourier transform of such exponentials. Graphs for 
a < 1 are in Ref A pp.305-3 11. 



Since the Fourier transform of the above correlation function from the FDT is a spectral density x 2 (co) for an extensive variable 
x corresponding to the concrete retardation, the relevant Levy distribution density is proportional to x 2 (co) d®, 



/(co) da> ~ x 2 (co) c/co. (2) 



This density is fractal at high frequencies, 



spectral density x 2 (co) — CO 1 a , 

susceptibility a"(co) ~ cox 2 (co) ~ co~ a , (3) 
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with the same exponent a; a" (go) is the loss part of the susceptibility (dynamic compliance) corresponding to x, CO the frequency; 
log CO is called the mobility. The general concept of fractality is restricted here to a power law expressing selfsimilarity. 

A Levy distribution with exponent a < 1 has neither a finite variance D nor a finite expectation E, 

D(co) = Var(co) -> °° for < a < 2, 

E(co) ~ J cox 2 (co) c/co^oo for 0<a< 1. (4) 

G)>0 

The second equation (0| leads to a preponderant component in the Levy sum for a Levy distribution considered as a stable limit 
distribution. One term of the sum retains a finite influence of order the "tilt" (1 — a) on the sum, irrespective that the number of 
its components tends to infinity. 

Are all these properties a direct consequence of the experimental Kohlrausch function (Q~|), especially for short times? 22 . Are 
additional arguments necessary that cannot simply be concluded from the existence of a Kohlrausch function or from present 
experiments at short times t or, correspondingly, at large frequencies CO? 

As indicated above, we try here an answer via a defect pattern from free volume V'(co), making higher mobility (logco) at 
places where more local free room is available for molecular motion. We shall introduce a Levy sum via separated partial 
systems or thermodynamic subsystems and try to explain the property a < 1 by some kind of instability. 

It is a general difficulty in understanding liquid dynamics by an approach that represents the complicated motion of molecules 
with "strong interaction" between many of them by means of statistical independence, being one of the basic assumptions to get 
a Levy distribution. 

The aim of this paper is a quasi thermodynamic approach to the dynamics of classical liquids starting from a robust 
mathematical and physical background for a fluctuating pattern with Levy defects. The mathematics used is based on the 
classical theorems of Darling and Feller. The mathematical development along their original proofs is accompanied by physical 
arguments, because the handling of a probability distribution with nonexistent expectation and variance (Eq. (0) is strange for 
physicists which are educated and experienced by only distributions that have normal (Gauss) limit distributions with finite, 
existing variances and expectations, in particular with no preponderant components in the limit sum. 



II. TERMINOLOGY FOR OUR LEVY SITUATION 

Aim. A glossarial introduction of some verbal concepts seems useful for the "thermodynamic" application of the Levy 
distribution in different fields. Calling f(x)dx = p(x)dx the density of a stable Levy distribution, then examples for the x 
coordinate of the sample space are: x = CO for the frequency of liquid dynamics, x — k for the wave vector of cosmic density 
fluctuation, and x = money for the economics of a Levy society ( 9 , p. 334-336). The discussion with the aid of Levy distribution 
is more suitable for the macroscopic, phenomenological thermodynamics (e.g. the state) and for susceptibilities near the 
equilibrium (modulus m(co), compliance j(co)); the discussion with a Levy flight 14 is more suitable for a molecular approach, 
e.g. the molecular diffusion leading to equilibrium values of a chemical potential. Warning: This Section is not a summary of 
the difficult probability theory of stable distributions^. 

Concepts. (1) Levy sum? is defined by 



S„ = Xi + ... +x„ 



(5) 
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where the Levy sum components X, are equivalent (cf. Point (2), below) independent random variables. The sum is suitable for 
the discussion of additive situations, e.g. for extensive thermodynamic variables, compliances, or for money in economics. (For 
the index i cf. Point (4), below). 



(2) Equivalence of the Levy components means that their distributions Fi(x) differ only by location ((3,-, b) and scale or norming 
parameters (at, a), e.g. for i= 1,2, 

F 2 (jc)=Fi(auc + Pi) ,ai >0. (6) 

Mathematically it is said that they are "of the same type". Another notion for two distributions U and V is: U = a\ + b. Often, 
mathematically, it is said that "...there are constants a, and (3; so that...". In the applications we need robust reasons for that: 
physical, cosmological, economical.... The "same type" is used for general arguments, the "parameters" are used for individual 
properties, e.g. for subsystems 

(3) Exclusiveness. Levy scaling. The former concept means that there are mathematically stable limit distributions (« — > °°) 
only for Gauss (a = 2) and Levy (0 < a < 2), with a the Levy exponent (cf. Eq. [7] and Point (4), below). For variations of 
components (small «), there are domains of attraction to the limit only for this two cases. It is important for the applications, that 
there is a broad spectrum of possibilities which tend to only the two sharply defined limits: Gauss or Levy, for the free volume or 
dynamics of liquids, or for the invisible hand of economics. We have then, for the n — > °° limit distributions, only one behavior 

of norming constants. In the form of S„ = c„X„ +y„, we get c„ = n 1 / 01 (Levy scaling). A useful concept is the scaled Levy sum, 

S„=S n /n l ' a . (7) 

This sum is existent for n — ► °°; the components X; are then equivalent to S„. A norming factor n- l l a is for "damping" the n ' 01 
"increase" of the original nonscaled Levy sum S„.This damping is much stronger than for Gauss (1 /^/n for a = 2). Examples 
are the sharp hierarchy of the ordered Levy sum for a < 1, Ref. 4 , or the non-Arrhenius behavior of the dynamic glass transition 
at low temperature (Section VI. D), cf. also Theorem 1, below. 

(4) Pluralism. This concept describes the great variety of the Levy distribution for applications. Two Examples. A. According 
to Eqs. ©, the Levy exponent a distinguishes three general cases (apart from some special cases for a — 2 and a = 1): (a) 
Expectation E and Variance D both finite: Gauss, a = 2; (b) E finite ("existent"), but D infinite ("nonexistent") for 1 < a < 2; 
and (c) both E and D infinite for < a < 1 ; 1 — a>0is called a "tilt". The (b) region is called hierarchy, the (c) region is 
called dictatorship (because we find a preponderant Levy component there) on the top of a sharper hierarchy. Example for (b): 
stock market with usually a « 1.5. Further. The Gauss case may be called democracy. The general equivalence of the Levy 
sum components is called primary democracy. — B. The index i of the components X, in Eq. (f5]l must be completed by further 
aspects: For the dynamics of liquids by temperature, pressure, kind of response after different disturbances; for the Levy society 
by the field (economics, politics, culture), and so on. The pluralism is a wide-ranging concept: Even for the same situation or 
the same defect, different response can have different exponents with different dependence on temperature, pressure, and so on. 
Such "individual" a exponents are not fixed by general arguments (cf. Section VI. C, below). 

(5) Fractality characterizes here the power laws of Eq. © and their selfsimilarity. Defining^ F (x) or F x (x) as the distribution, 

X °° 

f(x) or f x (x) as the density with F x (x) = J f x (Qa%— f(x) is continuous for Levy distributions — and 1 — F x (x) = J f x (£,)d% 

for large x coordinates as tail, the fractility according to Eq. (01 holds for the large x, i.e. for the tail: f(x) ~ x _1_a , with a 
the Levy exponent. Complementarity, for the characteristic function (e.g. for the simpler Laplace transform instead of Fourier, 

<p(A.) or <P/(A.) = / e~^ x f{x)dx = E(e _ )), we get a stretched exponential according to Eq. ([TJ with t = time for Fourier and 

o 

x = co frequency in a spectral density p(co) = /(co). The stretched exponential is more suitable for large times (or for less money 
in economy), where fractality is more suitable for high frequency (or many money in economy). The input of a spatial aspect 
(Eq. ( fT4l i. below) leads to a spatial unity, called "defect", with, for a < 1, a preponderant component in the midst, a fractal 
(hierarchical) center, and a stretched, cooperative periphery, all with the same exponent a. The size of the defect may be defined 
by a cooperativity N. 
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(6) Levy instability. Mathematically, the term "stable" means that, whatever the deviations in a finite sum of independent 
random variables in the domain of attraction are, the limit sum S„ stabilizes itself, i.e. S„/n l / a converges to a definite limit 
distribution for n — > °°. Exclusiveness for a < 2 leads then to a Levy distribution. In the applications, for illustrative or intuitive 
purposes, fractality for large x coordinate values, esp. with a preponderant component for a < 1, can be associated with a very 
special kind of instability, called "Levy instability" 9 . Example for liquids: a local breakdown to high frequencies due to a center 
of lower particle density, i.e. a local concentration of free volume. In the applications, of course, the high frequencies are 
very large but finite. "Very large" means a comparison of the center with the periphery, and £2 — > °° means a tendency, e.g. 
D (co) ~ f2 2 ~ a and E (co) ~ i2'~ a , formally with CO = Q. — > °°. The liquid example is described in Section III. A. 

(7) Levy situation is defined as an applicative situation (e.g. in physics, cosmology, economics) with a robust equivalence of 
the Levy components in the Levy sum and with robust reasons for a Levy instability for explanation of exponents a < 2, in 
particular of a < 1 when a preponderant component is to be expected. 

Examples. A distribution is symmetrical, if f(x) — f(—x); its expectation is zero, E(X) = 0. The spectral densities for 
classical (non-quantum mechanical) liquids are symmetrical, x 2 (— co) = jc 2 (co). Positive variables X > (i.e. the variables 
x > 0) are often physically motivated. For our (and many general) considerations about Levy distribution, X > is no serious 
restriction, because centering (selection of y„ in Point (3) above) is relatively free in the frame of equivalence. [To find e.g. the 
general expression for Levy distributions, depending on a and y (Ref. 3 , p. 580), no centering procedure is required for a < 1, 
while for a > 1 the natural centering to zero expectation suffices. In Particular, the influence of preponderant component is 
independent of y values. For construction of a concrete moleculare structure of a defect, however, centering may be important^.] 

For the sum of two independent variables, Xi +X2, the joint probability density f g {x\,»i) — fi(xi) ■ fiixi)- A common 
variable may be defined by the sum x=x\ + X2- The common distribution F(x) for P(Z, < x) is obtained by integration over the 
x space region x < x\ +X2 with measure dx\dx2- We get then for the sum density the convolution 



For the ratio of independent variables, X1/X2, we have e.g. a joint probability density f g (xi,X2) ~ *2 fi(x\X2) fi{xi) for 
F(x) = P{j* < x}21. The integration is more complicated than for the sum, because the ratio of the variables x\ jx2 < x must be 
used. The result is 



with the possibility of a partial rescaling procedure that affects f\ and f% differently: the product x\X2 vs. X2 alone. 

Be f\ a symmetrical Gauss density (a = 2) with variance Di = a\, and f% such a density with D2 = a|, then, for independent 
variables, the sum gives again a symmetrical Gauss density with variance D = a 2 = o\ + a 2 , (Gauss + Gauss = Gauss). Then- 
ratio, however, gives a Cauchy distribution, 




(8) 




(9) 



O1O2 



(10) 



K(G 2 + G2X 2 ) 



This is a symmetrical Levy distribution with exponent a = 1, i.e. with no existing variance (D = °°), and applied for positive 
variables x > 0, also with no expectation (E = 00); (Gauss/Gauss = Cauchy). The ratio favours large variances, but a Gauss/Gauss 
ratio cannot generate preponderant components. The latter are restricted to Levy exponents a < 1. 
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III. FLUCTUATING PATTERN WITH LEVY DEFECTS 



This section is to describe pictures for Levy instability and for equivalence of Levy sum components, i.e. a Levy situation for 
the fluctuating pattern of molecular liquids. 



To get a Levy instability, we consider a local breakthrough of molecular mobility (logco) as response to a virtual local 
concentration of free room for molecular movements. The aim of the concept of free volume (V') for the pattern is that a local 
concentration of free volume, i.e. a low local particle density, pushes the mobility there to high frequencies, in the limit to CO of 
order Q. 10 12 rad/s as for free rattling of particles in a sufficiently large molecular cage of next neighbors. Small local density 
reduction of order some percent is sufficient to rise co from e.g. 10 6 rad/s in the main transition to cage rattling, because the 
repulsive inter (and relevant intra) molecular potentials are steep. Locality of the breakthrough follows from a given free volume 
of the total sample in equilibrium. Lower local density here demands higher local density somewhere else. Lower density 
everywhere would finally lead to negative pressure. As the degree of divergence in expectation E and variance D is larger for 
smaller Levy exponents a, we may say that this exponent controls the Levy instability. 

To get a large number of equivalent small systems as components of a Levy sum we partition the volume of a subsystem into 
small parts, called partial systems: The large numbers needed for a robust dynamic partition for the limit distribution of spectral 
density x 2 (co) are thought to be numbers of random attempts for rearranging the molecular situation in the nanometer range 
there. We have many high-frequency attempts for a relevant low-frequency (co) rearranging in the dispersion zone considered; 
the number of attempts (as used for the limit) in the time interval 1/co is much larger than the number of affected particles. We 
think about molecular collisions which may cause the emission of a quantum /zco that can be catched by Nyquist's 25 transmission 
lines in a model for thermodynamic response, cf. also Ref.— , p. 270. 

Be careful to distinguish the two possibilities. The large number n for the spatial systems leads to a Gauss distribution for 
free volume with weak contrast in the pattern, the large number n for the attempts leads to a Levy distribution with exponent 
a < 1 for the spectral density with strong contrast. The relation between the two is explained in Step 2 of the Proof for the 
Representativeness Theorem in Section VI. D, below. 

We now assume that the frequency of slow rearranging (co) is locally controlled by the local free volume (V')- For given 
independent partial systems (index z) we assume, therefore, the existence of functions C0, (Vi) for each partial system, 



(no sum convention; Fig. 1). This assumption is justified by the high attempt rate. If the partial system is not too small, i.e. if 
the attempt rate inside is still much larger than the rearranging frequency, 



then we find a limit distribution for each partial system. Robustness of equivalence and partition for the domain of attraction 
and exclusiveness of Levy distribution allows to find a local control ( fTTT i and to identify the partial systems = 1, . ..,«,} with 
independent components of a Levy sum for free volume, 



A. Free volume. Partial systems. Levy instability 



C0i = C0/(^') 



(11) 



CO < attempt rate < Q. w 10 rad/s, 



(12) 



s'„ = ^'(coi) + ...+y>„). 



(13) 



B. Cooperativity 



A problem for intuition is whether the statistical independence of partial systems with infinitive expectation due to 
high attempt rates is sufficient for explanation of a phenomenon that is usually described by cooperativity or cooperative 
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FIG. 1: Local control of mobility (logco) by free volume (V 1 ) of partial systems, or by free volume of representative subsystems (index i, Eq. 

CD). 



rearrangement 26 . In the limit, each component has an equivalent dynamic distribution x 2 (co) [that is of the same type as all 
others, including the sum]. This is valid also for partial systems in the periphery ( TBl i of the defect. Their contribution to the sum 
(general susceptibility (O) may be small, nevertheless the movements in these partial systems may become very lively because 
of the infinite expectation (0J, E — > ». Any partial system can profit from the possibilities. It can e.g. take the advantage to use 
the high frequencies for an effect on the rearrangement, e.g. to become preponderant and get the center: All particles may take 
part in diffusion. 



Competition for cooperativity can so be simulated by independent attempts of partial systems. Expressed e.g. by the Gauss 
distribution for free volume, large fluctuation of one partial system i must be compensated by small fluctuation of some others, 
because the fluctuation of the sum becomes relatively small (1 j \fn). We get the chance to determine the extent of cooperativity 
(N mt ) from thermodynamic variables (Section VILA, below). For the Levy distribution with exponent a < 1 for dynamics, the 
competition of the interchangeable, equivalent Levy sum components leads to a preponderant component, probably connected 
with the diffusion step through the cage door, necessary for equilibrium in liquids. 



How can cooperativity be distinguished from collectivity, sometimes also used for the glass transition? Two types of 
collectivity are used for glass transition: (I) A common treatment of the high frequency c and a processes 20,27 (Fig. 3, below), 
and (II) A certain structure of thermodynamic phases (order parameter, clustering 15 ^). 



Let us here discuss type (I) collectivity. The separation of c and a process is related to two neighbored processes that both 
are at high frequencies and have direct mechanical relations. A Gotze ansatz can, therefore, use one formula for both, e.g. 
m(t) — vi cj)(f ) + V2§ 2 (t), where m{t) is a memory and (j)(f) is the relevant correlator, the mode of the mode coupling theory 
MCT This leads to a cusp bifurcation into a and c and to several scaling properties around the bifurcation, with reasonable 
experimental confirmation. Collectivity is therefore a suitable concept for the a and c processes in the relaxation chart of 
dynamic glass transition. 



What, however, about the crossover region from the a to the p and a processes (Fig. 3, below)? This is no bifurcation (Section 
VI. E below), which was assumed by MCT at the beginning. The rare experimental information about the crossove r 28 ' 29 shows 
two scenarios different from a bifurcation. It is discussed below (Section VI. E) that the defect pattern alone suffices to explain 
the existence and location of the Johari Goldstein p process. We see that collectivity is not needed for understanding the a, p, 
and (j) processes of the dynamic glass transition. 
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Remark. The above "inside" argumentation via partial systems for obtaining a Levy distribution is consistent to an "outside" 
proof via minimal subsystems (Fig. 4, Representativeness theorem Section VI. D). 



C. Johari's islands of mobility 

The concept of islands of mobilities is introduced by Johari 10 . In our context, we argue as follows. 

Since the free volume (V') has finite variance and expectation, there is no preponderant component in its distribution: we find 
a Gauss distribution for V'. For a dynamic spectral density x 2 (co) of the "same" partition in partial systems we obtain, however, a 
preponderant component due to the Levy instability (Step 2 of Section VI. C). A spatial concentration of free volume is, however, 
possible, if the interchangeable free partial or sub systems can spatially be redistributed. Some shallow "defect" in a free volume 
pattern (as basis for a mobility or frequency pattern from Eq. (fTTT i) can then be constructed by means of a general length/time or 
length/mobility scaling 11 for relevant modes (Fig. 2a and 2b), 




space < > space 

< > 

^ mt 



(c) (d) 

FIG. 2: (a + b). Defect pattern. If the white regions have more free volume, and the hatched regions less free volume, then (a) = islands of 
immobility that are not consistent with the general length/time scaling H41 , but (b) = islands of mobility are consistent 30 (dynamic hetero- 
geneity), (c). Sketch of the spatial defect structure. Calling the spatially separable contributions to the susceptibility oc"(co) ~ cox 2 (co) the 
"responsibility" [with j 2 (co) the spectral density via the FDT], we obtain from the above scaling and the Levy distribution a sharpened spatial 
picture for one defect. The one Levy distribution corresponds to the three parts of the defect: the influence of the preponderant component 
(for a < 1), the fractality of the center, and the cooperativity of the periphery, (d). The size of the defect in the pattern is determined by the 
average distance (Section VLB): The statistical independence of minimal representative subsystems for the mt dispersion zone determines a 
characteristic length h, mt via von Laue thermodynamics (Section VILA). 



defect periphery : large mode length ~ low mobility 
defect center : small mode length ~ large mobility 




This means that the co(V') function of ( fTTT i increases very sharply at the upper boundary for local free volume fluctuations. We 
get the picture of Fig. 2c. In the center of the defect we find the partial systems with much free volume, in the periphery those 
with less free volume. Transformation to the frequency (fTTT > means that the preponderant component related to the spectral 
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density is near or in the midst of the center for this island of mobility. In summary, we find more free volume in the center of 
the defect. 



Let us underline again, that the density differences across the defect are small, since the repulsive intermolecular potentials 
are steep. We get a dynamic heterogeneity with large mobility differences (but no or a small "structural" heterogeneity): 
a fluctuating pattern with temporary mobility defects. In terms of point (4) of Section II, the Johari Goldstein p process 
(Section VI. E) contributes to the dynamic dictatorship of the preponderant component. 



D. Description of the defect 

1. The defect embraces three properties of the Levy distribution (Fig. 2.c): (i) Fractality properties (O-© that care for the 
influence of large frequencies in the spectral density, (ii). A stretched exponential (fl} for the correlation function, (iii). The 
preponderant component of the Levy sum for spectral density. 

Assuming the above space/time scaling ( fT4b . then the Levy distribution stabilizes the Levy instability by connecting the 
center of the defect (fractal divergencies at high frequencies, preponderant component in the midst of the center, collection 
of partial systems with more free volume) with its cooperative periphery (long-time tail of relaxation from the stretched 
exponential, collection of partial systems with less free volume). Embracing by one distribution with one exponent a results 
in the correspondence between center and periphery (without any freezing-in). The statistical independence and equivalence 
of the Levy sum components inside the defect describes spatial aspects of the competition for cooperativity, e.g. for being the 
preponderant component in the midst of the center or for belonging to the hierarchy around the midst. 

2. We have three aspects of the preponderant component in the defect, (a) Induction of an extraordinary process for volume 
and entropy fluctuation, (b) The n — > °o limiting process for the Levy sum pushes the preponderant component in the midst of 
the center (Sections V.B-C). (c). The molecular picture is the diffusion step through the cage door and the promotion of the 
Johari Goldstein process. 

3. The size of the defect results from the length/time scaling ( TT4"i >. We find an opposite behavior of responsibilities (Fig. 
2c) on the way from a given defect to a neighbor defect. The responsibility in the periphery of the first defect decreases and 
becomes smaller than the increasing responsibility in the periphery of the second defect (Fig. 2d). We get an average distance 
in the pattern. The absolute length (in nanometers) can be determined by thermodynamics, if there are reasons to consider 
equivalent defects as some minimal representative units (Sections VI. C and VILA, below). 

The three points 1.-3. for the defects allow some visualization for the shaping power of Levy statistics in form of a fluctuating 
mobility or responsibility pattern in liquids. 



IV. RELATIONSHIP BETWEEN KOHLRAUSCH FUNCTION, LEVY SUM, AND LEVY INSTABILITY. THEOREM 1 

A Kohlrausch correlation function (Q~|i can also, in principle, be imagined without a relationship to Levy distributions, e.g. as 
obtained from a dynamic differential equation or even from a Hamiltonian. As mentioned in Section II, additional things are 
necessary to relate Eq. (Q]l with a Levy distribution for spectral density (x — co). This is expressed by the "if" in the following 
(RefAp.448). 

Theorem 1. For fixed Levy exponent a, < a < 1, the function (p( a ' (A.) = exp(— X a ) is the Laplace transform of a distribution 
f( w )(x) with the following properties: 
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(1) F^ a \x) is a Levy distribution; more precisely, if Xi , ...,X„ are independent random variables with distribution F a (x), then 
the normalized (scaled) Levy sum (Xi + .-.+Xn)//! 1 /" has again the distribution (x). 

(2) Fractality is obtained that may be normed by 

x a [l - F^(x)} 1 /r(l - a),x -► oo. (15) 

Comment. According to our basic assumption (Section I) we think here about spatially separable partial systems when 
components of a Levy sum are considered. For a physical Levy situation we need, for Theorem 1, a robust partition in 
equivalent partial systems, and, for Levy exponent a < 1, a Levy instability. In general, more trivially, it is the sum that allows 
exclusiveness of the Levy distribution to become a limit distribution for large numbers n. 



Proof of Theorem 1. Let us first of all recall that the probability density for a sum of independent variables may be 
transformed to a convolution (0 of the component densities that are (Fourier or) Laplace transformed in a product. Therefore, 
in the limit we get for any stable distribution an infinitely divisible distribution. The scaling property n l l a of stability (1) follows 
from infinite divisibility of the Kohlrausch function (stretched exponential (fl]i): (cp^ a ')"(5t) = ^ a \n l / a 'k).— The function cp' a ' 
is completely monotone. Since <p( a '(0) = 1, the measure F^ a '(x) with Laplace transform <p( a ' has a total mass 1. 



The proof for part (2) is obtained along the following line. A positive function L defined for < x < °° varies slowly at °°, if 
for every fixed x 



L(ax) 

-V^ -> 1 for fl^oo. (16) 
L(a) 

Such functions may be used to describe the domain of attraction of a probability distribution. The limit of Eq. ( TToT l is only 
fulfilled by power functions (the basis for our fractility) which is, after a special selection of the probability by the Gamma 
function of the tilt (1 — oc), T(l — a), expressed by ( fT31 >. End of the Proof. 



V. PHYSICAL UNDERSTANDING OF THE PREPONDERANT COMPONENT FOR LEVY DEFECTS 



Let us recall Feller's (Ref. 3 , p. 172) mathematical heuristic. "Consider, for example, a stable distribution ... with a < 1.... 
The average (X\ + ... +X„)/n has the same distribution as Xin~ 1+1 / a , and the last factor tends to °°. Roughly speaking we can 
say that the average of n variables is likely to be considerably larger than any given component X^.. This is possible only if the 
maximal term M„ = max[Xi , ...,X„] is likely to grow exceedingly large and to receive a preponderating influence on the sum 
S„". Section V is firstly to bring some analysis used and is secondly to describe the preponderant component in Levy defects 
more precisely than in Section III. 



A. Darling Lemma 



To calculate the influence of a preponderant component as expressed by the Levy exponent a < 1 we need a method for 
handling the maximal component (M„) in a Levy sum (S„). This is demonstrated along the original Proof for following 
(Darling's 2 ) 

Lemma. The Laplace transform cp^(A-) = ^„(A.) of the random variable Z„ = S„/M„ (which ratio should characterize the 
above influence) is 



l 

UV=ne- x j (P J e~ Xy f (yp) dy)"' 1 /(P) dfi 





(17) 
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with (p, y) dummy variables for the probability density / of equivalent components {X,}. 

Proof. We start from the equivalence of the Levy components. There is no loss in generality by assuming Xi = M„, since 
each X,- has a probability of l/« of being the largest term (and P{xj = xj} = for i ^ j since the distribution is presumed to be 
continuous). The joint density f g (index g) is then 

f n/(Pi)/(P2).-/(P„) 
/ g (p 1 ,p 2 ,...,p„) = <^ ifPi=max{Pi} (18) 

I otherwise 

Then the relevant expectation ^„ = 



o 

x 

/ l ... /' e -H72+-+Yn)/P >< 





x/(Y 2 ).../(Yn)-/(P)-^Y2...^P 

•J-l 




-^/P/(Y)^Y> /(P)dP- (19) 







The Lemma follows from rescaling of the dummy variables: Pa!y d(y$). The upper limit 1 in the inner integral comes from 
y< p, i.e. from the maximality of the component for p. Note the different variables, p and yp, for the densities / in the Lemma 
dl7l . typical for ratios of random variables (|9). 



B. Influence of the preponderant component. Theorem 2 

This subsection is to calculate the expectation of the random ratio: Levy sum over its maximal component, Z„ = S„/M„ with 
variable (observable) for n ^ °°. 

Theorem 2 (Darling 2 ). The expectation of Z„ for the limit n — > °° with a Levy exponent < a < 1 is 

E(S„/M„)^l/(l-a). (20) 

The Proof is divided in two parts: The limit is first explicitly calculated for the limit distribution itself, and in a Comment the 
participation of the domain of attraction is discussed. We put for the inner integral of the Laplace transform (p^"' (A.) = ^„ (A.) = 
E(e~ xz ") from the Lemma, Eq. (fTTT l, 

l 

V (P) d ^p| e -^/(YP)^- (2D 
o 

Note the different variables in exponent and density. In the limit (n — > °°), only the tails of the Levy distribution are important 
(cf. the Comment). Hence we put for their fractility 

l-F(x) = c/x a (22) 

with a constant c = c(a, y) where y is ascribed to centering. There are no further parameters beyond (oc,y) in the general Levy 
density. Eq. ( 1221 is a consequence of Eqs. (Q3) and (1161 1. Then, from arithmetic manipulation in direction of the tails, such as 
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we get 



(1-F(7P))-(1-F(P)), (23) 



1 

¥ (P) = 1 - (1 - - xj e-** [(1 - F(tP)) - (1 - F(P))] dy. (24) 



The integral yields 



itegral d ^ f -X | <T»*[...] rfy = |a j {e^-l. 
o o 



c/y 



Defining the number 



we obtain 



Rescaling now with the variable 



o 



The wanted expectation value follows as 



(25) 



(26) 



¥ (P) = l-(l-ch) C 7P a (27) 



v = c/$ a , (28) 
(which eliminates c and therefore any centering y in c), gives in our limit 

m = — i — — • (29) 



E($(X)) = -^'(0) = -^/^(X = 0) = 

d , ^ _ ct . 1 

= -lim — (l-X + X =- . (30) 

x^o dX a - 1 1 - a 

The influence of the maximal term on the limit Levy sum is calculated as 1/(1 — a) in the ratio S„/M„. The only remaining 
parameter is the Levy exponent a. 



We may say, that the influence of the maximal term, the preponderant component, is of order (1 — a), (the tilt in cosmological 
terms). Compared with Gauss, where the influence of the maximal term remains infinitely small, we get e.g. for a = 1/2 
the surprising result, that the maximal, the preponderant component has approximately the same influence as the rest of the 
infinitely many components. 



Comment. The domain of attraction is handled by the deviation from Eq. ( fToT l, using the symbol o(l) which means of the 
order less than 1, 
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We pursue this symbol along the proof by the chain due to rescalings, 

o(l)->o(l-F(P))->o(l/fi). (32) 
From ( f2Tb we get for 1 — F(y$) after yp — > y rescaling 

1 -F(7p) = (1/7°) (1 -F( Y )(1 (33) 

The integral d25t yields then the 



integral = (l-F(P)) a J (e~**- 1) y-^dy+ofl -F(fi)). 



(34) 



The shift into the tail is therefore 



¥ (p) = 1 - (1 -00(1 -F(P)) -F(p)) (35) 

with the same (j)i as before, Eq. (f26b . Since |\|/(P)| < 1 for bounded (3, the early portion of the integral for £,„(k) is negligible for 
sufficiently large n. This can be seen after introduction of a new variable instead of d28l l, now 

v = n(l-F(P)). (36) 

Hence 

n 

&,(>•) [(l--(l-h) + vo(-)) n - 1 dv, (37) 

J n n 
o 

and we get again 1/(1 — a) in the limit. End of the Proof for Theorem 2. 



We see from the Comment chain (1321 that the preponderant component really pushes the play into the tail's fractality 
asymptote. The large number of participants (attempts in our model) pushes the preponderant component into the midst of the 
defect center; Fig. 2c is a picture for the shaping power of the Levy statistics. 

Remark. Absence of a preponderant component would mean S„/M„ — > °°. This is prevented by the o(l/n) pushing into the 
defect center and the finiteness of the remaining integral due to the fractality y _1_a in the denominator of d34l l. This leads to the 
finiteness of —£'(o) in ( f30l >. 



C. Verbal understanding of a preponderant component in a defect 

This brings some comments to the Points (4) and (6) of Section II and continues the defect discussion of Section III. The Levy 
distribution is characterized by a large number (n) of statistically independent equivalent random components in a sum (inside 
example: the partial systems of a subsystem; outside example: the subsystems of a larger subsystem) and a Levy exponent a 
that describes some Levy instability, larger for smaller a. Let us remark, that the equivalence of the Levy sum components 
means that they are interchangeable near equilibrium. They can be interchanged between the periphery, the hierarchy, and the 
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preponderance of the center. 

In a spatial picture of a a < 1 defect (Fig. 2c above), the preponderant component is in the midst of the fractal center 
of the defect. This component is pushed in the center by Eq. d32l . The defect center collects the components with more 
influence around the preponderant component: the sharp "hierarchy" below the "dictator". The other components are in the 
cooperative periphery where the components with temporarily less influence on the sum are collected, although all components 
are equivalent by having infinite expectation (E — > °o). The Levy distribution stabilizes the defect by connecting the fast 
fractal center with the slow cooperative periphery; there is a correspondence between center and periphery. The cooperativity 
is stimulated (encouraged) by the equivalence in expectation. For infinite expectation, the competition of the equivalent and 
independent components is stimulated by the property that a considerable part of influence, of order the tilt (1 — a), on the total 
sum is realized by only one, the preponderant component, interchangeable with the others in the equilibrium. 

In other words, the "actual events" are pushed in the center with the preponderant component. This are therefore fast actions. 
In liquids, the preponderant component in dynamics may induce the diffusion step through the molecular cage door as basic 
process to maintain the equilibrium. In the midst of the center there is only action for one, the preponderant component. [The 
"dictator" can fill the midst because of his infinite "personal" expectation, E — ► °°.] This is mathematically complemented by 
the absence of (this) one component outside, which absence makes the integrals of Eq. < fTvT > suitable. 



Let us conclude with an example for nonequilibrium. If we have an underlying trend, in liquids e.g. a cooling rate 
f = dT jdt < 0, then we may ask, which parts of the defect firstly cannot longer follow. Due to the general length/time scaling 
( TBI relevant for the spatial structure of the defect, it is the slow periphery that first falls out of equilibrium. The trend leads 
unavoidably to incrustation, to a frozen vault structure (Ref.— , p. 386), while the center continues a time with responsibility to 
outer dynamic contacts. 



VI. DYNAMIC GLASS TRANSITION IN MODERATE MOLECULAR LIQUIDS 

A. Dispersion zones in a mobility -(reciprocal) temperature diagram 

The term "dynamic glass transition" is synonymously used for the dynamics in liquids. For classical molecular liquids of 
moderate complexity ("moderate liquids") we find a surprisingly general relaxation chart, i.e. a log CO — 1/T diagram (Fig. 3). 
There is a characteristic arrangement of dispersion zones (zones of higher dissipation) with a typical width of one or a few 
frequency decades, if crystallization can be prevented. The symbols and names of these zones are listed in the figure caption. 
This picture can be discussed as a consequence 9 of Levy defects for the main transition (mt = a plus a process) and for the 
Fischer modes ((])). 

The concentration of mobility (logco) into zones can be explained by the defects. The width of the imaginary part of the 
susceptibility (dynamic compliance a"(co)) is, because of a"(co) ~ COx 2 (co), determined by the embracing Levy distribution for 
x 2 (co). The corresponding a"(co) graphs for a < 1 show (Ref. 9 , p. 307 and p. 31 1) an approximate half width of one decade in 
frequency co, divided by the Levy exponent a, 

Alogco = (1.05 ±0.05) /a. (38) 

[This equation is valid independently on the spatial dimension of the defect model. The width of the a" (logco) peak depends 
only from the Levy exponent a.] Large mobility distances between the zones in the relaxation chart care for large attempt rates 
in partial systems as used for the above inside treatment. 

The non-Arrhenius behavior of the main transition mt is, in principle, connected with the increase of its cooperativity N„u at 
low temperature (Fig. 6, below) via a Levy scaling equation 



d\og®/dT ~N}„ 



(39) 
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FIG. 3: Relaxation chart for moderate liquids; log CO the mobility, T the temperature. Dispersion zones: a the high-temperature process, a 
the cooperative process, (mt = a plus a = the main transition or main process = The dynamic glass transition), (3 the local Johari Goldstein 
process, C the a — a crossover region, and (]) the Fischer modes. [Boson peak b and cage rattling c are not further discussed.] The conventional 
glass temperature T g corresponds approximately to the intersection of the a process with the experimental 10 millihertz isochron. T c is the 
crossover temperature. More complicated substances such as polymers or liquid crystals may have additional dispersion zones. 



where N mt is the number of molecular units in the mt defect. This equation is related to the existence of a bulk modulus mt being 
the "reciprocal" of an invariant free volume v/ ~ 1 /my, 

Vf=dV'/dlog<i>, (40) 
invariant against the partition of a larger system into partial systems (Section VI. D). 

A Levy treatment of partial systems inside one minimal subsystem (Ref.— , p. 327, see also Section VI. D, below) is explicitly 
based on our main assumption to consider spatially separable components of a Levy sum. Further development 30 leads to an 
extreme smallness of mt for small Levy exponents a at low temperature, mt (a < 0.4) = o(mf(u » 1)) m.f(T = T c ). The 
large curvature of the mt in the relaxation chart follows from ( |39l l with a thermodynamic equation along the mt, dV' ~ dT: 
d\og(x)/dT ~ d\og(a/dV' ~ 1/vf ~m/ increases more than Arrhenius, d\og(a/dT ~ 1/T 2 . The equation 

mf-vf-N- 1 '" (41) 

is called control equation. This actual operational freedom m/ becomes narrow inside large units, irrespective of the dictatorial 
control by a preponderant component. 



Remark. The extreme smallness of the operational freedom ;«/ below a < 0.4 implies an exhaustion of the mt process at 
low temperature where, in the equilibrium, not enough free volume can be organized by increasing cooperativity of the defect 
periphery. The mt process is then successfully "attacked" from distributed "accidental" loss effects (Ref. 9 , p. 178), and the 
control by equation (|4TT > is expected to be lost at a surprisingly sharp 30 exhaustion transition. 



B. Gedankenexperiment for finding a minimal subsystem 



If external large heat reservoirs and the imagination of walls between subsystems are given up and if, for dynamics, the 
subsystems are defined by statistical independence in a given dispersion zone, then the larger subsystems are free in partition, 
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and all subsystems are freely fluctuating in all relevant thermodynamic variables. We obtain freely fluctuating subsystems. 

The following gedankenexperiment is for the definition of representativeness. Consider a sufficiently large freely fluctuating 
subsystem having a certain spectral density x 2 (co) (e.g. for fluctuating (free) volume in the main transition and the following 
increments for thermodynamic variables across the mt). Divide this subsystem into two halves. If they also are sufficiently large, 
then both have the same non-extensive (e.g. density-like) variables. It is said that they are representative. Divide again, till the 
subsystems become too small for representativeness. We find finally a minimal representative freely fluctuating subsystem for 
the given dispersion zone: shortly a minimal subsystem. 

Representative and minimal means that this subsystem contains one defect and has a size of order the defect size (Fig. 2d). 
The defect labels the minimal subsystem as a thermodynamic unit. The unit has three corresponding parts: preponderant 
component, center, and periphery, united by a Levy distribution. [Since the only true unity of the defect is the one preponderant 
component, we may say that this component makes the minimal subsystem to a thermodynamic unit: to "the molecule" of mt 
fluctuation of the liquid.] 

Representative subsystems have the same relevant thermodynamic variables as the macroscopic systems, with all its 
fluctuations. This means that also the minimal subsystem has a temperature T that can completely be defined via Carnot and 
Kelvin, and has an entropy S that can completely be defined via Clausius; both T and S fluctuate, Sr 7^ 0, 85 7^ 0, also for the 
minimal subsystem when it is in the nanometer range. The systems from the gedankenexperiment plus arbitrary combinations 
of them are expected to be candidates for equivalent and robust components in an appropriate Levy sum for thermodynamics of 
liquids. 



We conclude from the above gedankenexperiment that the shape of spectral density x 2 (co) in all representative subsystems is 
the same. We expect additionally from the equivalence of all partitions into a sum of representative subsystems that for additive, 
i.e. extensive variables the shape of x 2 (co) is fixed to be a density of a Levy distribution. We show, so to speak, that the bare 
existence of minimal subsystems in liquids implies the defect. As an example, we think firstly on x 2 (co) = AV' 2 (co) as a density 
for free volume (V 1 ) fluctuations in the main transition (mt). 

Representativeness theorem. (Ref. 9 , p. 237). Consider the classical (non-quantum mechanical = symmetrical) spectral 
density for stationary fluctuations of free volume, or a corresponding extensive caloric variable, in the slower part of the mt 
dispersion zone (logco (rad/s) < 1 1) of a representative freely fluctuating subsystem. This spectral density is a Levy distribution 
density f(x)dx with frequency measure dx = do) and a Levy exponent a < 1. The exponent a depends on temperature T, 
pressure p, substance, and kind of response. 

Remark. The appropriate variables for measurement are susceptibilities, e.g. dynamic ones as function of frequencies 
(OC*(co) according to Eq. (01; i.e. moduli ra(co) for intensive and compliances j'(co) for extensive thermodynamic variables). The 
thermodynamic increments for e.g. the extensive (free) volume follow from integrals across the dispersion zone (DZ), 



Defining the volume compliance by B = — dV /dp, [B] = m 3 /Pa, we get for the FDT (the "measuring equation" 9 p. 269 ff.,— ) 



C. Representativeness Theorem. Pluralism 




(42) 



DZ 



AV 2 



k B TB, 



(43) 



with a correlation function AV 2 



(t) in the time domain, 
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AV 2 (t)=k B T (fl(f)-flequii), (44) 
or with a spectral density AV 2 (co) in the frequency domain, 

AV 2 (co) = k B TB"{(S)) /too, (45) 
with B* (<b) = B'(co) — iB" ((£>), real part minus imaginary part. 

Proof. The proof is divided in five steps. 

1. Let us recall some terminological foundations of probability theory. Consider the random variable Xi = V[. The 
observable, Vi > or vi > for V[, means a coordinate in the sample space Br, i.e. a variable for a "conceptual experiment" of 
probability. [For thermodynamics, this experiment is described by the FDT]. The distribution Fi(vi) for the random variable 
is defined by the probability P to observe a value Vi < Vi, Fi(vi) = P[Vi < Vi]; the coiTesponding density is denoted by /(v), 

b 

P[a < v < b] = J fi (v\)dv\. Analogous definitions are used for the other subsystems, {V^ ■■■,Vj, ...V,,'}. 

a 

For the sum of independent free volumes V{ of two subsystems, S' 2 = V[ + V^, we have ^52(52) with a convolution (O for the 
density, 

fsi(s 2 )= jMs2-v)f 2 (v)dv = f 1 *f2(s 2 ). (46) 
The variance and expectation for any free volume exist, e.g. 



E(s 2 ) = J S2fs2(s 2 )ds 2 . (47) 

The distribution for S' n as function of the variable s„ tends, therefore, to a Gauss limit distribution. In general, we have a Gauss 
bell curve for /y(v) = /(v). 

2. Interested in dynamics of liquids, we must change the variables from volume to frequency, v — > v(co). Consider the 
representative inverse function co(v) of Fig. 1, e.g. via the inside treatment of the local breakthrough of mobility, our Levy 
instability for a < 2. Practically, substituting co = Q. w 10 12 rad/s by co = 00, we find a continuous function diverging at some 
finite, maximally accessible free volume inside the minimal subsystem, < °o: 

K>(v)-^°o for v = v max . (48) 
The expectation value of frequency may therefore diverge, 



E((0)=E((Q(v)) = J C0(v)/(v)rfv->oo, (49) 
and we expect a new distribution, fy (co) = /(co), different from Gauss. 



One may ask, whether the general shape of the new, dynamic density /(co) depends on details of the co(v) function. The 
answer is no, since the convolution property (|46T > for many components v; pushes the distribution for their sum into the domain 
of attraction of a limit distribution. Then from a divergence (|49l we expect a Levy distribution with Levy exponent a < 1, 
whatever the details of the continuous co(v) function beyond the property (|49l are. The details of co(v) can, however, influence 
the particular value of a. It is interesting that changing the measure may imply a change of the exponent a (cf. the Corollary). 
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3. What is the appropriate (physical) density function /(co) for the stationary fluctuation of an extensive variable? This is 
the spectral density, /(co) = x 2 (co) (0, because the corresponding variable measures are equivalent as indicated by the integral 
Eq. ©. 

J x 2 (v)du = J f(v)dv = AV,' nn (50) 

where AV„' jr is the (free) volume step increment for crossing the main transition (mt). The "additivity" of jc 2 (co) is also ensured 
by its equivalence to the compliance dV /dp or dV /dT via the FDT 

4. From (l46i l for free volume we get, after the substitution f(v) — > /(co) = jc 2 (co) for each component of the Levy sum, 

X?(<o) = (X 2 *Xl)(<6). (51) 
Hence, for general «, we have for any sum of minimal subsystems 

X 2 n ((0)=X 2 *...*X 2 ((0). (52) 

Physically, the general frequency variable co follows from the co-identity of the FDT, measuring a general frequency for any set 
of representative subsystems and the same frequency CO for any activity, 

{co,} - co. (53) 
Since d52l is a sufficient condition for a Levy distribution, we get such a distribution forx 2 (co) with measure c/co. 

In the Nyquist model 9 25 for the FDT (mentioned in Section III. A), the co-identity is related to the indistinguishability of the 
/zco quantums in the transmission lines between "sample" and "apparatus": different activities correspond to different selections 
of quantums from these lines. 

5. A thermodynamic argument for exponent a < 1 follows from the periphery of the defect. The Levy distribution connects 
the center of the defect (fractality) with the periphery (tail of stretched exponential, Fig. 2c). The long-time retardation is in the 
periphery. The simplest case for retardation is a "thermodynamic" or "quasi-stationary" one (in the terminology of Landau and 
Lifshitz, Ref. 32 , §118) 

dx/dt ~ —x , x — > for t — > oo. (54) 
This gives a Debye (exponential) decay corresponding to a Cauchy distribution (O (a = 1) in the frequency domain. 



The cooperativity of the periphery induces a spectrum of such thermodynamic decays. This means stretching the retardation. 
For Levy distribution we get necessarily the Kohlrausch function (HJ, i.e. a Levy exponent a < 1, holding for the total defect 
because of the correspondence between center and periphery. This ends the Proof, if the question of variables wherefrom the 
exponent a may depend (and physically depends) is delegated to the 

Corollary. The only free parameter in the correlation function for fluctuations of extensive variables in the main transition 
of moderate classical liquids is the Levy exponent a < 1 . The exponent a is representative, but remains non-fixed for minimal 
subsystems in the one-nanometer range. It is therefore expected to be influenced by the structure and the responsibility (Fig. 2c) 
of the molecules and of the molecular structure of the Levy defect. This includes the dependence on the thermodynamic state 
(temperature T, pressure p, composition x,...). We have got 



a = a(activity, substance, T,p,x, ...). 



(55) 
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Activity means the different kind of response for extensive variables, e.g. volume, entropy, dielectric polarizability, incoherent 
or coherent dynamic neutron scattering, and so on. 

Six Remarks. (1) The Representativeness Theorem shows that in moderate liquids the outside approach to the Levy 
distribution via representative subsystems is consistent with the inside approach via partial systems. 

(2) The Gauss distribution for free volume from Eq. ( f39l > excludes the picture, that there is a partial system with a preponderant 
contribution of the free volume itself. But this Gauss distribution does not exclude a spatial concentration of such partial systems 
that do have more free volume than others. The general space/time scaling (Fig. 2) can generate a defect with more free volume 
from these partial systems in the center. 

(3) Does the correspondence between center and periphery mean that there is no diffusion step through the cage door in a 
Levy liquid without a preponderant component? 

(4) Long range correlation as from lattices in solid state physics and long range quantum correlations in liquids may prevent 
a reasonable partition into smaller and smaller independent subsystems, so that representativeness does not longer remain a 
reasonable concept: No spatial Levy defects there. 

(5) In other words, the Corollary expresses some kind of pluralism of different activities. Although the different activities 
are collected by the molecules, in a way, in the same dispersion zone, and although always the same molecules of the minimal 
subsystem are involved, the activities are different in shape and location across the zone. This is related not only to the value of 
the Levy exponent a, but also to modifications from molecular particularities in the center (e.g. proximity of NMR signals to 
mechanics) or in the periphery (e.g. nonlocal properties of shear response). Each activity for itself is related to Levy, but each 
has different Levy exponents and possible modifications, even in the same (p,T) state. The shape of the limit distribution has 
different "faces". Although frequently confirmed by experiments, such a pluralism of a theoretical approach seems difficult to 
be accepted by the glass transition community. 

(6) A drastic example for the pluralism is the above variable-change from free volume v' to frequency, co(v'). For/(v) = f(v') 
we get a Gauss distribution (a = 2), and for the spectral density /(co) = x 2 ((o) we get a Levy distribution with a < 1. Such 
findings may have consequences on the statistical physics in liquids (Section VII, von Laue approach). The problem is hidden, 
because Gauss behaviour (a = 2) may be expected for all thermodynamic potentials that may be found in the nominator of 
Gibbs exponents, from similar arguments as for free volume. 



D. Derivation of invariant free volume and control equation 

The Levy treatment of Section VIA is based on an inside treatment of cooperativity, inside one minimal subsystem i.e. inside 
one defect (Section VLB). Consider independent partial systems with a fixed average volume. The size of the subsystem is then 
proportional to the number m of such partial systems, N mt ~ m. Consider the frequency CO,-, i = 1 , m, as a transition probability. 
Statistical independence is then expressed by a common product, 

(£> g = const x coi x ©2 x ... x co m . (56) 

From ( fTTT i and Fig. 1, we have 

CD, = C0, (V/) , same i ; (57) 

from additivity of free volume V 1 we have 
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V' (subsystem) = V[ + + ... +V' m . (58) 
These three equations define a functional equation yielding the general solution 

v! + Vi + + V 

co(subsystem) = <D M exp(— - — — -) (59) 

v f 

with com a prefactor depending of the partition and vj a scaling parameter of the minimal subsystem not depending on the 
partition. This vy is called the invariant free volume. The cooperativity is now described by exchange of free volume between 
the independent partial systems balanced by the relatively small fractional volume fluctuation 5V/V of the total minimal 
subsystem ("minimal coupling"). 

Intuitively, the decrease of operational freedom m/ (Section VI. A) is connected with smaller free volume. Comparing directly 

the scaled Levy sum £X,//? 1//a (Eq. (0, cf. also Theorem 1, above) with the exponent in Eq. (1591 , we would get Vf ~ N„, a , 
i.e. Vf would increase with the cooperativity. This means that vy cannot be directly proportional to this freedom. To analyse 
the situation (free volumes V/ also decrease with lower temperatures) we are led to susceptibilities near equilibrium, being the 
direct measures from the FDT (Section VI. C). We have then dynamic compliances y'(co) that become larger for softer materials 
(high temperature, small cooperativity AO, and moduli m(co) that become larger for harder materials (low temperature, large 
cooperativity). That is, the operational freedom, decreasing for low temperatures, increases with the modulus. Using the 
rheological equation j'(co) ■ ra(co) = 1, and putting j ~ Vf ~ N^J t a , we get for the operational freedom m/ ~ m(co) ~ 1/j ~ vT 1 , 

i.e. really Eq. (|4TT >. [In the book 9 and Ref 30 this relation was used without any derivation.] Since the general Levy scaling by 
n 1 /™ is for Gauss consistent with the fluctuation variance (— , p. 172), the scaled Levy sum seems to be the only possibility to get 
a general damping of the operational freedom with increasing cooperativity N. Reducing with Gauss (a = 2, nif ~ N~ l l 2 ), we 

get m r — m//m/(Gauss) ~ 2 " 2a ) e g TO; ^ N m ^ 2 for a = 1 /2, i.e. a considerably stronger restriction of freedom than 
from Gauss. 



E. Johari Goldstein (3 process 

Consider falling temperatures. The p process "emerges" in the a-a process crossover region C of the relaxation chart Fig. 3. 
According to the idealized version of the mode coupling theory MCT 20 , the cage door of next molecular neighbors (for the a 
process), necessary for the diffusion in liquids, is closed at the critical temperature T c of this MCT. In our Levy defect approach, 
T c is substituted with the crossover temperature: Because of the generality of defects, esp. some continuity of molecular motion 
in the cooperative periphery (see Fig. 6, below), Johari's islands of mobility (Section III B,C) survive the crossover with a 
modification 33 . The increasing cooperativity of the a process, located in the periphery of the defect, leaves the cage door open 
for diffusion. The position of the crossover region in the relaxation chart is then determined by crossing of the a process with 
a virtual molecular-mechanical process for (3, finding free room enough in the center of the defect. From this point of view, 
the parameter change for the mt at the crossover 33 means a change from a defect without (T > T c ) to a defect with the inside 
mechanical process (T < T c ). [The parameters of the a process and the mechanical core of the p process must suffice some 
relation to get a crossing. There are glass formers where no crossing is observed, i.e. with "no p process"—]. 

A p process is approximately of Arrhenius type and is accompanied by a "continuous" background from the periphery. 
The total formal variance of the p process can be affected (i.e. broadened) by a distribution of molecular activation-barrier 
properties. A concrete p activity can be rather small, if the mechanical process does not generate the corresponding fluctuations. 

The main issue is the relation of the p process to the a process below the crossover region. Wherefrom does the free volume 
concentration necessary for the above activation process come? Although the free volume of the bulk becomes extremely small 
at low temperatures, so that the a process needs more and more cooperativity for perpetuating the diffusion steps through the 
cage door? A short answer is the phase continuity of the liquid state so that cooperativity must increase to ensure the fluidity of 
the liquid, even at the cost of long times (large viscosity) at low temperature, cf. Eq. ( |39l ). The longer answer from our approach 
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is that free volume enough can be collected in the center of Levy defects of the cooperative a process (with the preponderant 
component for dynamics in the midst). The correspondence between center and periphery then guarantees that increasing 
cooperativity of the periphery for falling temperatures (Fig. 6, below) ensures the continuous existence of the mechanical 
process by the preponderant component for a < 1, also when the total free volume becomes shorter. [Let us repeat: In the 
idealizing MCT closing of the cage door is possible; in our approach the periphery cannot end its action to keep the cage door 
open.] The alternative 35 , that the p process is a precursor or local step of the a process which can occur anywhere 36 , cannot 
answer the above wherefrom question. In our defect picture, the anywhere would only be right in the sense that the Levy defect 
can occur anywhere in the pattern, but anywhere outside the defect center would not be probable. 



F. Fischer speckles 



One of the greatest surprise in glass transition research was the discovery, that the Debye Bueche 37 inhomogeneities in 
frozen glasses correspond to a slow dispersion zone (]) in the relaxation chart Fig. 3. These Fischer modes have e.g. times 
about 10 7 times longer and lengths about 50 times larger than those of the main transition process (mf)i&. It follows from the 
Representativeness Theorem that there must be also defects in the fluctuating pattern for the (j) process zone. These defects are 
to be identified with the experimental 17 Fischer speckles. 



The volumes of partial systems for the inside approach to the Glarum defects of the main transition are smaller than the 
volumes of the mt minimal subsystems, whereas for the Fischer speckles the partition into just such subsystems (as "partial 
systems" for (j)) is used (Fig. 4). 



fast dispersion zone: 
main transition 
(mt process) 



P 



center of 

Glarum Levy defect 




slow dispersion zone: 
()> process 



f 



) speckle 



FIG. 4: Fischer speckles ((]) speckles). Two Levy situations in moderate liquids. Left for the main transition mt, right for Fischer modes (]). The 
small squares are the minimal subsystems for the main transition. 



The Fischer modes (j) are a direct consequence of the main transition mt. No structural thermodynamic assistance ("cluster") 
is necessary for the construction of the length and time ratios between (j) and mt processes (Fig. 5). These ratios can be estimated 
from alone the two-defect pattern of Fig. 4. We have to look for a physical mt process in the liquid that eats its way through 
large space X and time t ranges: diffusion process. How probes the fast mt process the slow (j) process speckles? What about the 
reliability of the short-range mt process for the long-range (j) process? 



The mt process defects are rather dense, their average distance is of order 1 nm (see Fig. 6, below). The molecular diffusion 
is therefore step by step, from one Glarum Levy defect to the next one along the "diffusion line". We expect locally, at any 
step, a modification of the "normal diffusion" by the Levy fractality at the cage doors: "Levy diffusion". This defines the local 
log t j log A, slope, and an "analytical continuation", step by step along the line gives the required reliability of the defect behavior 
over the full distance between the speckles. The Representativeness Theorem includes also the large subsystems between the 
Fischer speckles. — Conversely, the fast mt process should also be "probed" by components of the slow (j) process. We have got 
two aspects: 



Seen from the fast mt process. For instance, dynamic neutron scattering • shows that the a process can be described by a 
sublinear diffusion whose exponent may be related to the Kohlrausch (Levy) exponent a of a certain a process relaxation. In the 
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process (|> process 

with tilt (1-a) 

FIG. 5: Tilt construction for Fischer <|) process length (t,§) and time (t^) in a diffusion plot, log? vs log A, 2 , using the relevant tilt (1 — a) of the 
Glarum Levy defects of the main transition (mt). 



frame of pluralism (Section VI. C), this exponent may be different from a dielectric exponent (cf. Eq. (l55ll) 40 . The corresponding 
Levy diffusion is described by a slope from 



co — > co a , t A 2 /" (60) 

Seen from the slow (]) process. Consider the flaring up (1) and fading away (0) of the (j) speckles at different space and time 
"points" as some kind of (1,0) spin diffusion from virtual, spatially fixed spin carriers. Since different minimal subsystems for 
the mt process are independent, this spin diffusion is normal: 



r-A 2 . (61) 

We estimate the required (j) /mt ratios from the triangle in Fig. 5. We find the nearest neighbors of (]) process defects (speckles) 
at the upper (j) corner of the triangle, since the Levy diffusion, to be effective, must be faster than the normal diffusion. The basis 
e is the log-time difference 8 between normal and Levy diffusion at the mt process length scale, X « t, mt . The Levy diffusion is 
at some time t in the short time tail (fractality) of the mt process spectral density, i.e. in the molecular cage door being the center 
of the Glarum Levy defect. The normal diffusion is near the average time of the mt process, t av : 



e = log(t av /t) mt . (62) 

From Fig. 5 we see that the time and the length ratios ((|)/mf) are mainly determined by the tilt (1 — a). Defining T„„ = t av , we 
get the slope difference as ( 1 — a) /a and 



i ^ a i T <e a 

log ir- = wr, r£ , log — = e- (63) 

c> 2(1 -a) v, 1-a 

For small tilts, 1 — a <C 1, we obtain the slope difference directly as the tilt (1 — a), and log^/c^,,,,) w e/2(l — a) as well as 
log(T^/x mr ) — a). Large values for the (j)/mf ratios can easily be obtained, also in the experimental range for moderate 
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molecular liquids (Ref. 9 , p. 349). A minimal subsystem for the process as claimed by the Fischer speckles contains a large 
number of Glarum Levy defects. 



VII. VON LAUE THERMODYNAMICS FOR LIQUIDS? 
A. Characteristic length of the main transition 

The characteristic length is a concept for an estimation of the defect size from thermodynamic bulk measurements (Fig. 2(d), 
Sections III.D and VI. C). This length is decisive for an Experimentum Crucis of the defect pattern. 

The usually applied molecular Gibbs statistics rests on a virtual transfer of the temperature T from a macroscopic heat 
reservoir, transferred without temperature fluctuation to the subsystems, 87" = 0. This transfer does not seem reasonable, if one 
tries to eliminate external experiments (the conscious observer of the Copenhagen program) from thermodynamics or quantum 
theory of independent subsystems. 

Statistical independence of Landau subsystems is rested on negligible potential energy between the foreign molecules 
according to the Boltzmann factor exp(— E pot (p, q) /ksT) in the Gibbs probability. This does not seem appropriate for dynamics 
that is based on large attempt rates allowing much smaller independent subsystems. 

It seems more reliable to substantiate a thermodynamics of liquids by time or mobility integration of dynamic compliances 
and dynamic moduli from the fluctuation dissipation theorem (FDT) that can be derivated from quantum mechanics without 
external experiments (RefJi, p. 227ff j 25 ' 31 , cf. the Remark in Sec. VI. C). This leads, for (T,S; p,V) as relevant thermodynamic 
variables, to the von Laue thermodynamics including temperature fluctuation &T also for the above minimal subsystems in the 
one-nanometer range (Refs i 32 > 41 §112, and 42 Chap. 14). 

A thermodynamics based on the dynamic compliances (oc"(co) ~ COx 2 (co))— without any restriction of fluctuation of 
thermodynamic variables, similar to the von Laue thermodynamics — may possibly be associated with a statistical basis 
different from the more "static" Gibbs distribution. In particular, the separation into mechanics (E(p,q)) and thermodynamics 
(ksT) of the Gibbs distribution may be lost. Possibly, there exists a molecular, microscopic von Laue distribution with different 
consequences near the defects, related to the Levy distribution as representative of minimal subsystems and reflecting the 
pluralisms (different a for different activities) more explicitly. 

From the von Laue fluctuation formulas we get directly 43 a volume that can be interpreted as volume of the minimal subsystem 
when applied to the main transition mt. We start from one defect as the smallest representative thermodynamic unit (Section 
VI. C). In the general von Laue fluctuation formula for temperature, AT 2 = kgT 2 /Cy, the heat capacity Cy is an extensive variable 
[J/K]. Having N mt molecular units in our minimal subsystem, we get 

iV m ,(vonLaue)=^^k ) (64) 

where now the small cy letter is the specific heat capacity [J/kg-K]: A(l /cy) mt is the step of l/cy across the mt, R is the molar 
gas constant, and Mq is the molecular mass of what is considered as molecular unit. The characteristic length is obtained from 
the average volume needed by the N mt molecular units; N mt was called "the cooperativity" for the defect of the mt. 

We discuss now the extensivity (being the expression of the Levy sum). Put v the number of minimal subsystems in the 
sample; i.e. Eq.d64b is for V = 1. To get the extensivity of Cy, we multiply the heat capacity of one minimal subsystem with 
V and get then for the temperature fluctuation of the sample Ar 2 ~ 1/v, AT ~ l/v^ as expected. But how to get (8r) 2 r , 
the temperature fluctuation of the minimal subsystem, from a bulk measurement? The Representativeness Theorem makes the 
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relevant responsibility (oc"(co) ~ coAx 2 (co)) indepent from the system size, i.e. the width Alogco of the mt dispersion zone (DZ) 
of Fig. 3 does not become smaller for large subsystems or samples. We get, across the DZ, a temperature difference 



5r » (dT/d\og(x)U ongDT • 5logco. (65) 

To interprete this ST as hT mt , we must assume that for partial systems of the minimal subsystem, the distributions for 
temperature and relevant mobility fluctuations are equivalent. This is a fiuctuative interpretation of Boltzmann's macroscopic 
temperature-time equivalence used by rheologists for a long time. [This equivalence is an approximate symmetry of susceptibil- 
ity a = a(logco,r) curves across the mt with respect to a logco <-> T exchange, called also rheological simplicity. Interpreting 
in this context log co as a thermodynamic variable equivalent to a temperature T, this symmetry can be thought to be inversely 
constructed from fiuctuative increments for both, 8logco and 5T. This confirms our above assumption, in particular for relevant 
susceptibilities a, e.g. for temperature modulus or for entropy compliance ~ dynamic heat capacity, C((0,T). In the FDT 
as thermodynamic equation for an experiment, this equivalence is mediated by the Planck radiation formula in Nyquist's 
transmission lines, cf. Fig. 3.6b of Ref. 9 ; a detailed preliminary discussion of this issue is also in 9 , pp. 263-268. A critique of 
this argumentation is in — .] 



Eq. ( 1641 is from the fluctuation of an intensive variable, the corresponding susceptibility is a modulus. From a Gibbs fluctua- 
tion of energy, (AE 2 ^ V = ksT 2 Cy (with the restriction 8r = 0), we get from an analogous treatmen t 45 ! 46 

RT 2 

AW (Gibbs) = - . (66) 

Mo (Acy) m oTgjnt 

This is from the fluctuation of an extensive variable, the susceptibility is a compliance; hT gjnt is the transformation interval after 
correction for partial freezing around the thermal glass temperature T g . The ratio is unexpectedly large, of order 



N mt (Gibbs) / c x 



N mt (Laue) \Ac 



(67) 



where Ac is the step height across the mt of the heat capacity used for the estimation 47 , and bT 2 w bT 2 mt which assumption is 
not a point of issue. For e.g. Ac/c = 10 per cent we get the ratio of 100 (!) 



The two experimental cooperativities as obtained from dynamic calorimetry between mHz...kHz frequencies (3(0 method) 
are compared in Fig. 6 48,49 ' 50 . We see three things. 



(i). The Gibbs cooperativities are much larger than the von Laue cooperativities. The related Gibbs lengths correspond about 
to the diameter of the cage from next molecular neighbors plus parts of the second shell, as expected for Landau subsystems 
from the energetically motivated independence due to the Boltzmann factor. 



(ii). There is no indication for any singularity (cusp or so) at the crossover temperature ' 9 (for the substances where T c is 
included) that could be expected from a critical temperature from closing the cage of the pristine = idealizing model for mode 
coupling theory MCT-22i2I. 



(iii). The von Laue cooperativities from calorimetry increase continuously with falling temperature. The increase corresponds 
to a continuous decrease of bulk free volume. Some of the cooperativities are, for high temperature above the crossover, 
T > T c , smaller than one molecular cage from the next neighbors and suit, therefore, to the extraordinary entropy fluctuation 
accompanying the diffusion step through the cage door induced by the preponderant component in the midst of the defect 
center. Sometimes, the cooperativities are of order N mt ~ 1: one molecular diameter near T c . Some bend in the crossover 
region is also observed. For low temperature, the increasing cooperativities seem to indicate the increase of the defect periphery 
to collect enough free volume for the center to maintain the cage door open also for short free volume, as required by the 
Representativeness Theorem. 
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FIG. 6: Comparison of von Laue and Gibbs cooperativities (Eqs. J 64b and ( 1661 )) of the main transition mt as a function of reduced temperature 
x (T c the crossover temperature, Tq the Vogel temperature from a VFT (Vogel Fulcher Tammann) extrapolation to log CO — » — °o; x= 1 means 
T — T c ). The full experimental data are in Refi^, see also Refsi 4 V° and the references cited therein. The colors are: blue for poly (n alkyl) 
methacrylates, green for other polymers, and red for small-molecule substances. For polymers, the particles are the monomeric (repeat) units. 
The numbers are: 2, 3, 7, 8, 10 for poly (n alkyl) methacrylates: 2 - n propyl, 3 - n butyl, 7 - n pentyl, 8 - n hexyl, 10 - n decyl; 4 - random 
copolymer poly (n butyl methacrylate-.staf-styrene) with 2 percent styrene; 13 - polystyrene, 16 polycarbonate, 19 diglycid ether of bisphenol, 
20-poly [(phenyl glycidil ether) - co - formaldehyde]; 21 - benzoin isobutyl ether, 24 ortho cresyl glycidyl ether, and 25 - a fulven derivative 
(TP-CPBO). The bars mean an error estimation of the single absolute values. 



The experimental characteristic lengths from the von Laue approach are, therefore, considered as size of the Glarum Levy 
defects. 



B. Experimentum Crucis 

Is there an Experimentum Crucis that can decide between using von Laue or Gibbs thermodynamics and, if von Laue is the 
winner, would indirectly prove the Levy defects in the main-transition (mt) dispersion zone of the dynamic glass transition 
(Fig. 3; see also Ref.S p.187, 50 )? Since the local density effects expected from the cage doors for diffusion are small due to 
the steep repulsion potentials of the particles, it seems at present difficult to detect the defects directly in the structure factors 
used e.g. for the evaluation of dynamic neutron or photon scattering. Additionally, the inter/intra molecular potentials used for 
adjustments are also uncertain. 



Therefore, an indirect experimental proof via characteristic lengths is suggested. Dynamic scattering of relevant contributions 
from mt process allows to construct a raster 51 of wave vectors Q, mobilities logcoo, and isotherms T = const. Transferring 
the maximum frequency C0q C for the heat capacity peak of Cp((0, T) at different temperatures from dynamic calorimetry (DC) 
into the logcoo-^ lines of this raster, the resulting length Z, = 2k /Q from scattering can be compared with the characteristic 
lengths from the calorimetric cooperativities, Eqs. d66b or d64i i. Since their ratio: Gibbs/von Laue is large, in particular for high 
temperatures (Fig. 6), even a rough comparison can give a decision. If we obtain von Laue lengths, then this result is considered 
as a certain experimental proof of the defects (see below). 
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What about the physical relevance of calorimetry for neutron scattering in the frame of pluralism (different exponents a)? 
Roughly, by thermodynamics, scattering corresponds to some density fluctuation, i.e. to a volume, and calorimetry to an entropy 
fluctuation. The von Laue fluctuations are related by thermal expansion: ASAV = kgT (dV /dT)p. This defines a certain rela- 
tionship between them. The general basis for comparison, however, is the co-identity of the FDT ( 9 , p. 279). Take thus dynamic 
heat capacity and the diffusion part of dynamic scattering as the two responses that are to be compared. Then the comparison 
of maximal frequency coo of the two compliances Eq. Q seems reasonable, since the frequencies of different activities (even 
for different Levy exponents a = ockww) are identified by the co-identy. The Kohlrausch function from the relevant intermediate 
scattering function, S ml (Q,t) ~ exp{ — [f/ T KWw(2)] a } must be transferred to a dynamic compliance, a"(co). The maximum 
frequency coo can then be calculated from Ref. 9 , p. 306: 

ln(cooTKww) = 0.60607 (a -1), (68) 

where In means the natural logarithm, coo is the angular frequency in rad/s, with CO = 2nf, f the frequency in hertz. We get 
<*>o(2)- The reduction to frequencies avoids further comparison of times from different activities. 

Consider a fictitious example following partly Ref. 40 (Fig. 7). From dynamic neutron scattering (DNS) one can separate a 
relevant intermediate scattering function S (Q,t) that is informative about diffusive (Eq. [60b and related relaxation processes. 
This is one of our response. Calculation of co DNS (Q, T) from a Kohlrausch fit with Tkww as a function of wave vector Q allows 
to construct a raster of isotherms in an logC0Q NS vs \ogQ plot. This raster allows then to locate the maximum (logcoo, T) pairs 
from dynamic calorimetry (as the other response) in the diagram Fig. 7 (A) (by using co^ = co DNS from the co-identity of 
the FDT, possibly after some extrapolation of the DNS raster to lower temperatures). Then the abscissa of the A points gives 
the Q vector that is to be associated with dynamic calorimetry, from which the length \ = 2n/Q can be compared with the 
characteristic lengths from ( l64j i and (|66] l, (i^Jp ) 3 = N mt x average volume occupied by one particle. 




DNS 



FIG. 7: Fictitious diagram for the Experimentum Crucis. The maximum frequency values C0q (T) of dynamic heat capacity Cp(co, T) for 
different given temperatures, (A), [located by the raster of isotherms, and using co DC = co DNS from the FDT] in a log frequency (logCflg NS ) — 
log wave vector (log<2 DNS ) diagram for that part of dynamic neutron scattering which is informative about diffusive and corresponding 
relaxation processes. DC the dynamic calorimetry, DNS the dynamic neutron scattering. The log Q abscissa values of the DC triangles A give 
qualitative measures of the lengths (L, = 2n/Q) that can be associated with DC and can be compared with the calorimetric characteristic lengths 
from j64\ and l |66t . Fig. 6. Details see text. 

In Ref. 40 , DNS for polyvinylacetate PVAC is compared with the peak of e"(co), the dielectric loss. The authors get a length 
of about 1 nanometer (from Q = 0.65 A -1 ). In our experience with DC and pluralism, the trace of e" in the relaxation chart 
is usually located not too far from the C" trace. The extrapolation of von Laue cooperativities to x = 1 (crossover) gives the 
order of one or a few particles, i.e. the same length order. The Gibbs length would give cooperativities N,„ t of order 100 - 1000, 
corresponding to lengths of about 2-5 nanometer, i.e. (much) larger. This preliminary comparison actually indicates von Laue 
thermodynamics. 



27 



In a way, the defect model is directly supported by DNS for polyisopren o 38 ' 39 in a large Q interval (0.1 A _1 ...5.0 A -1 ) and 
time interval (Tkww = 0.03 ns...200 ns). The raster of Fig. 7 corresponds to a large "homogeneous" region (Levy diffusion) 
for small Q values (large lengths \ = 2k /Q), which is continued by a "heterogeneous" region for large Q values (smaller 
lengths). The Tkww times there are above the raster extrapolation. The crossover between the two regions is at Q w 1.3 A -1 , 
corresponding to i; as 5 A, close to the maximum of the static structure factor, S(Q). This observation can be interpreted 
by our defect model, Fig. 2c. The homogeneous region corresponds to the fractal center of the defect. The large lengths 
can be explained by the analytical continuation as used for the Levy diffusion to the Fischer modes (Section VI. F). The 
heterogeneous region corresponds to the diffusion step through the cage door. A responsibility part of this region corresponds 
to the preponderant component of the Levy sum, being, so to speak, an "inhomogeneous" contribution from the homogeneous 
Levy distribution. The other part is from the heterogeneous molecular contributions in the construction of the cage door up 
to large Q (i.e. down to small lengths). The observation of the homogeneous part supports one of our basis assumption for 
the Levy statistics, that the large numbers needed for limit distributions come from the many dynamic attempts in the partial 
systems, and is not directly the number of particles participating at the cage door. 

Distinguish two results from the Experimentum Crucis. (A). Decision between von Laue or Gibbs thermodynamics for 
minimal subsystems in liquids, and (B). The values of the bulk characteristic lengths per se, in nanometers, indicating the spatial 
extent of the entropy or temperature responsibility in the defects. 

It should be tried to widen the reliable experimental ranges of dynamic neutron scattering along the mt dispersion zone to 
lower frequencies than today really available (about several 10 7 Hz) and of dynamic calorimetry to higher frequencies (than 
the today 10 4 Hz), so that both methods can overlap for the same substances e.g. in the crossover region of the mt (usually in 
the 10 6 Hz range, cf. Ref, 9 , p. 222-223). This would decide which alternative for the characteristic length would be consistent 
with scattering. The typical glass transition crossover frequency can be reached by dynamic neutron scattering in the next 
years, if e.g. the impressive series of experiments by Dieter Richter's group (e.g i 38 i 39,52 ) is continued to lower frequencies. 
Reliable dynamic calorimetry in the megahertz range, however, remains to be done 53 . Dynamic photon scattering can probably 
be applied in the a-process dispersion zone at low frequencies, a few frequency decades above freezing-in around the glass 
temperature T g . Experimental data from dynamic calorimetry is available her o 47 ! 48 . 

The decision for von Laue would become a conclusive experimental argument for Levy defects, if our treatment from defects 
to von Laue thermodynamics has a robust reversal. This includes three main steps: (1) Freely fluctuating subsystems should 
give a sound thermodynamic basis for the components of a Levy sum with Levy exponent oc < 1 (Representativeness Theorem). 
(2). In particular, such equivalent subsystems should give a thermodynamic basis for the way from Kohlrausch correlation 
functions to the Levy sum (Theorem 1). (3). We need a sound explanation to have a Levy instability in the liquid giving a < 1 
(free volume picture Fig. 1 and Theorem 2). I think we can explicitly find the robust reversal in these steps. 

There are general consequences for thermodynamics, if the temperature fluctuation at so small, one-nanometer lengths could 
be confirmed by the Experimentum Crucis. Since we do not have a molecular spatio-temperature field in this length range, i.e. 
no T(r,t) there, we cannot illustrate temperature fluctuations by additional hydrodynamic arguments using forces and fluxes 4 - 2 -. 
Instead, 8r must be included in the many possibilities of the molecular origin for the random components of a Levy sum which 
is agitated by the high attempt rates of dynamics. Instead of hydrodynamics, the Levy defect forms some Levy diffusion. 
Temperature fluctuation is there a necessary part of the dynamic molecular play in liquids. 

Teaching thermodynamics, we could leave the thermodynamic limit for the definition of temperature. If 8T ^ would be 
confirmed, then the Gibbs distribution comes from a "static" restriction: we start from a more general thermodynamics for 
subsystems with bT ^ and come to the limit N — * °° with 8r = 0, e.g. by means of a heat reservoir. The use of this restriction 
as a basis of the general definition of temperature (as sometimes advocated 54 ) would not longer be reasonable. Temperature 
T is essentially different from an energy AE, irrespective of the frequent occurrence of an equation of type T = AE/ks- In 
liquids, such an equation comes from Planck's black body radiation of the quantums AE = h(ti in Nyquist's transmission lines 
for a thermodynamic measurement, /ico = kgT. In the FDT, this equation connects temperature fluctuations and frequency 
fluctuations for a specific response, as used above. 

In summary, an Experimentum Crucis which confirms the lengths calculated from von Laue thermodynamics would overcome 
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the nonfluctuation of temperature — i.e. brings a new (the old) definition of temperature in liquids beyond (and before) the 
Gibbs definition with the latter's heat reservoir without fluctuation. The consequence for liquids could be that mechanics and 
thermodynamics cannot be separated into denominator over nominator, (E(p,q) /ksT), as for Gibbs. The "statistical mechanics" 
could be influenced by a thermodynamics with Levy defects in the liquid. Then the observables, in particular the dynamic 
susceptibilities, cannot exactly be calculated from mechanical computer simulation plus Gibbs thermodynamics alone. 



VIII. CONCLUSIONS 

Most of the qualitative attributes for the arrangement and properties of the dispersion zones (a, a, p, (j)) in the relaxation chart 
of the dynamic glass transition (as a synonym of classical liquid dynamics) can be explained by a weak underlying fluctuating 
free-volume pattern with strong dynamic Levy defects. Eight Examples are: 

(1) The Levy defect explains the concentration of the main transition (mt) dynamics into a relatively narrow dispersion zone 
with a half-width of about (1 /a) decades, with a < 1 the Levy exponent. The spectral densities for extensive variables are Levy 
distributions, their correlation functions are stretched exponentials. 

(2) The spatial position of the Johari Goldstein p process in the pattern is the center of the Glarum Levy defect of the a 
process; the continuous part of the p process comes from the cooperative periphery. 

(3) The location of the crossover region in the relaxation chart is defined by virtual crossing of a quasi mechanical molecular 
barrier mechanisms (Arrhenius process underlying the p process) with the main mt process being not an Arrhenius process. 

(4) The emerging von Laue thermodynamics allows an experimental determination of the characteristic length for the main 
transition by bulk dynamic calorimetry. This length is the size of minimal representative subsystems and characterizes the size 
of the Glarum Levy defects (including its preponderant component, its fractal center and its cooperative periphery) and their 
mutual average distance in the pattern. 

(5) The preponderant component of the mt Levy sum for the defect may induce an extraordinary molecular event of dynamics: 
The diffusion step through the cage door opened in the cage of equivalent next molecular neighbors. 



(6) The non-Arrhenius behavior of the main transition mt (= a + a process) is connected with the increase of cooperativity 
N mt at low temperatures. 



(7) The position of the dispersion zone for the Fischer modes relative to the main transition in the relaxation chart is 
determined by the diffusion-relevant Levy exponent a from the mt Glarum Levy defects without the need to refer to any 
thermodynamic collective structure formation beyond the Levy distribution (no "clusters"). The Fischer speckles are the Levy 
defects of a Levy sum from a large number of minimal subsystems. 

(8) At low temperatures (or mobilities), the main transition is probably exhausted. The operational freedom in the periphery 
for large cooperativities N mt becomes so small that spatially distributed, random traces of foreign processes can take the control 
over from the main process. This restricts the Levy exponent from below, a(mt) > 0.4. 

The way from the widely observed stretched exponential (Kohlrausch function) to a Levy distribution for the spatial 
fluctuation pattern with preponderant components as midst of the defects is explained by a succession of theorems which 
are partly known before. Their original mathematical proofs are partly modified and completed by physical arguments and 
comments. Two things are additionally needed for this way: Mathematically, (1) a sum of Levy sum components and (2) an 
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explanation of small Levy exponents a < 1. This corresponds physically (ad 1) to a sum of spatially separable statistically 
independent subsystems or partial systems and (ad 2) to a Levy instability from a local breakthrough of mobility. The latter 
is characterized by more free volume in the center as well as a spectrum of relaxation in the cooperative periphery of the 
defect. This "inside explanation" is confirmed by an "outside" one from the Representativeness Theorem for thermodynamics 
of liquids, operating with a sum of representative minimal subsystems. An Experimentum Crucis for the whole argumentation 
is suggested: the comparison of characteristic lengths from thermodynamics (dynamic calorimetry) with lengths from dynamic 
neutron (or photon) scattering. 



The Levy distribution gets a "shaping power of statistics" if it is imbedded in certain additional circumstances. For the 
dynamics in classical liquids, e.g., the following is sufficient: (a). Many equivalent independent spatially separable partial 
systems for its establishment as a limit distribution, (b). An instability whose stabilization by the limit distribution is different 
from Gauss (exponent a < 2, in particular a < 1 for the existence of a preponderant component in the Levy sum), and (c). A 
general time (mobility)-length scaling of relevant fluctuation modes to get a fluctuating space pattern with defects. The latter 
are thus considered as the results of the shaping power. 
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